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Variational  Aj>proximationt  to  the  Diffraction 
by  Circtilar  and  Elliptical  Apertures 


by 

Chaang  Huang 

Cruft  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 


Abstract 


Different  methods  are  reviewed  briefly  for  attacking  the  problem  of 
the  diffraction  of  a plane  electromagnetic  wave  by  a circular  aperture  in 
a plane  screen  which  is  infinite stm^y  thin  and  perfectly  conductixig. 
Specifically,  the  variational  method  derived  by  Levine  and  Schwinger  has 
been  used  to  evsduate  the  transmission  'coefficient  of  circular  and  ellipti- 
cal apertures.  For  the  circiLiar  aperture,  a high-order  vector  trial  fvmc- 
tion  with  undetermined  frequency -dependent  coefficients  is  chosen.  By 
using  the  stationary  property  of  the  eaq>ression  for  the  transmission  coef- 
ficient, equations  for  the  undetermined  coefficients  are  derived.  These 
equations  are  solved  to  give  a first-order  approximation  of  the  transmis- 
sion coefficient,  and  the  numerical  values  are  compared  wilth  previous 
results.  For  elliptical  apertures,  a zeroth-order  approximation  of  the 
transmission  coefficient  is  evaluated  using  a one-component  trial  ftmction. 
Numerical  results  are  given  for  ellipses  with  nodnor-to-maior-axis  ratios 
of  1/2  and  1/3. 


I. 

Introduction 

The  problem  of  the  diffraction  of  electromagnetic  waves  by  a circular 
aperture  in  an  infinite  plane  conducting  screen  has  attracted  the  attention 
of  many  authors.  Their  work  is  reviewed  briefly  in  this  section. 

Moglich^  tried  to  work  out  the  exact  soltition  for  the  problem  of  dif- 
fraction  by  a ciz'cular  disk  (which  is  the  Babinet  complement  of  a circular 
aperture)  by  expanding  the  rectangiilar  components  of  the  Hertz  vector,  v, 
of  the  incident  and  scattered  fields  in  series  of  oblate  spheroidal  wave  func- 
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tions  and  then  imposing  the  boundary  condition  on  ir  at  the  surface  of  the 
disk  to  determine  the  unknown  coefficients  of  the  expansion  for  the  scat- 
tered Hertz  vector.  His  solution  was  incorrect  because  be  did  not  take 
into  consideration  the  singxilarity  at  the  rim  of  the  disk.  This  defect  was 
pointed  out  by  Meixner,  who  reconsidered  the  problem  by  expanding  the 
Debye  potentials  into  series  of  oblate  spheroidal  wave  functions,  and  later 

4 

improved  his  own  resiil^  with  the  collaboration  of  Andrejewski  by  em- 
ploying ^e  Hertz  vectors.  Their  solution  includes  the  complementary 
problenr,  in  the  Babinet  sense,  of  diffraction  by  a circular  aperture.  The 
analogous  exact  solution  to  the  more  general  problem  of  the  elliptical  disk 
awaits  some  future  study  of  ellipsoidal  wave  functions. 

Although  exact  solutions  for  an  arbitrarily  shaped  disk  or  aperture 
are  not  feasible  at  present,  a number  of  different  formulations  ave  avail- 
able for  obtaining  approximate  results  over  specific  frequency  ranges. 

The  application  of  these  methods  to  the  problem  of  a circular  disk  or  aper- 
ture permits  a comparison  with  the  results  of  the  exact  theory  and  gives 
some  idea  of  the  degree  of  accuracy  to  be  eiqiected  from  a particular 
approximation.  It  is  to  be  emphasized,  however,  that  the  significance  of 
these  methods  lies  in  their  applicability  to  disks  or  apertures  that  are  not 
circular  in  shape. 

To  the  extent  that  they  are  solutions  of  Maxwell's  equations  all  these 
formulaticns  of  the  problem  of  diffraction  are  equivalent.  When  explicit 
results  are  desired  for  a partievdar  case,  however,  the  necessary  approx- 
imations give  differing  degrees  of  accuracy,  depending  upon  the  frequency 
range  that  is  of  interest  and  the  formulation  that  is  used.  Approximations 
which  have  proved  tiaeful  in  the  past  can  be  classified  as  follows: 

1.  The  lowVfrequency  (static)  app roximation : The  aperture  or  disk 

is  taken  to  be  smali  compared  to  the  wavelength;  ka  = 2ira/X«l  (where 
a = characteristic  dimension  of  the  aperture  or  disk)  . 

a.  Retardation  is  neglected,  and  the  result  gives  the  first  few 
terms  in  an  expansion  in  powers  of  ka. 

b.  The  incident  field  vectors  are  taken  to  be  constant  over  the 
small  aperture  or  on  the  surface  of  a small  disk. 
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Using  these  approximations  Rayleigh,^  Tai,^  Bethe  * ^ and  Bouwkamp® 
have  studied  the  problem  of  diffraction  by  small  circvilar  disks  or  aper- 
tures, The  results  obtauned  by  these  investigators  are  accurate  only  for 
frequencies  which  satisfy  the  condition  ka«  1 . 

2.  The  high frequency  (optical)  approximation:  The  aperture  or 
disk  is  taken  to  be  large  compared  to  the  wavelength;  ka»l  . These 
approximations  are  of  the  Kirchhoff  type. 

a.  The  tangential  component  of  the  magnetic  field  vanishes  on 
the  dark  side  of  the  conducting  screen. 

b.  The  tangential  component  of  the  electric  field  in  the  aperture 
is  equal  to  the  incident  electric  field. 

Formulations  which  make  use  of  these  approximations  have  been  dis- 

10  11  **12 
cussed  by  Stratton  and  Chu,  Schelkunoff,  Silver  and  Ehrilich,  and 

Levine  and  Schwinger.  Among  these  four,  the  analysis  of  Levine  and 

Schwinger  has  the  advantage  of  a relative  analytic  simplicity,  and  is 

14 

recapitulated  in  a previous  report,  where  its  application  to  annxxlar, 
elliptical  and  rectangtdar  apertures  is  worked  out. 

13 

In  contrast  to  the  above  methods,  Levine  and  Schwinger  have  derived 
vector  integral  equations  using  dyadic  Green's  fimctions,  and  they  have 
shown  how  to  calculate  the  far-zone  diffracted  fields  and  the  transmission 
coefficients  of  apertures  in  terms  of  variational  principles  related  to  the 
integral  equations.  Thus,  by  choosing  an  appropriate  trialiunction  for  either 
the  tangential  component  of  the  electric  field  in  the  aperture  or  for  the.  mag- 
netic field  on  the  back  side  of  the  diffracting  screen,  an  accurate  result  can 
be  obtained  for  the  transmission  coefficient  without  making  any  further 
asstimptions.  In  this  way,  they  have  computed  the  transmission  coefficient 
of  a circular  aperture  for  normal  incidence,  using  two  different  trial  fimc- 
tions  for  the  tangential  component  of  the  electric  field  in  the  aperture. 


9 

*Bethe's  solution  has  been  modified  by  Bouwkamp. 

*‘*Silver  and  Ehrlich  have  solved  >4axwell's  equation  by  a Fourier  oper- 
ational methpd  and  obtained  formulae  xlnr  the  near -zone  fields. 

**‘*Copson^  ^ has  also  derived  the  Iniegral  equations  for  the  diffraction  of 
disks  and  apertures,  but  since  he  usee  scalar  Green's  function,  his  for- 
mtilatioh  is  very  cumbersome. 
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The  firit  of  theae  functions  waa  constructed  from  a aingle^component 
electric  field  with  no  angular  variation  and  with  a radial  dependence  cor> 
responding  to  that  of  the  static  solution.  It  results  in  a transmission  coef- 
ficient that  is  accurate  only  at  low  frequencies.  The  second  trial  function 
is  the  complete  static  aperture  field  obtained  by  Boxtwkamp.  It  gives  an 
improved  result  which  is  still  accurate  at  moderately  high  frequencies 
where  the  Bouwkamp  solution  itself  is  not  valid. 

Since  these  trial  functions  are  frequency -independent,  they  give  xinsat- 
isfactory  answers  in  the  middle  frequencies.  In  this  study,  a frequency- 
dependent  trial  function  will  be  constructed  for  the  circular  aperture,  and 
it  will  be  seen  that  a very  high  degree  of  accuracy  can  be  obtained  with 
only  moderate  complexity.  Finally,  the  variational  procedure  will  be 
applied  to  elliptical  apertures  using  a single -component  trial  function  anal- 
agous  to  the  first  one  used  by  Lievine  and  Schwinger  for  the  circular  aper- 
ture. 


n. 

The  Circular  Aperture 

The  Stationary  Expression  for  the  Transmission  Coefficient 

The  problem  to  be  considered  is  that  of  a plane  electromagnetic  wave 
incident  on  an  aperture  S^^  which  perforates  an  infinitesimally  thin,  per- 
fectly conducting  plane  screen,  S^.  A rectangular  coordinate  system 
( X,  y,  z)  is  chosen  so  that  S^  and  S^  lie  in  the  plane  z = 0 as  shown  in 
Fig.  1. 

The  plane  wave  is  incident  in  the  half-  space  z<0  and  is  described  by 

ginc^r)  = ^ expv(  ikn.  r ), 

Inc  /N  C ^ ) 

H “*^(r)  = hexp(ikn*r), 

where  r is  a position  vector  in  space,  fi  is  the  unit  propagation  vector,  ^ 
and  £ are  unit  polarization  vectors,  and  k is  the  wave  number.  In 
gauBsian  units  the  unit  vectors  are  related  by 

^ = ^x^  , and  s £.h  = 1 
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The  harmonic  time  dependence  e 
throughout  this  report. 


•i«^ 


is  omitted  here  and  will  be  omitted 


The  transmission  coefficient  of  an  aperture  is  defined  as  the  ratio  of 

the  transmitted  energy  flux  per  unit  aperture  area  to  incident  energy  flux 

per  unit  area.  A stationary  expression  for  this  quantity  has  been  derived 

13 

by  Levine  and  Schwinger;  it  is 

(h*y^ £x£  {p)exp(-ikn- p)dS)  (h*/^ zxE  (p)exp(ikn>  p)dS) 

\ ^ -n 


/ 


S X E^(p)-  r^°>(p  ,p')  S X E_^(p)  dS  dS' 


(2) 


where  p is  the  position  vector  in  the  aperture  and  A is  the  area  of  the 
aperture,  are  the  aperture  electric  fields  produced 

by  incident  plane  waves  with  propagation  vectors  n and  -n  respectively. 


both  having  the  magneto  polarization  vector  h 
space  dyadic  Green's  function  and  is  given  by 


r^®^(p,p')  is  the  free- 


(3) 


4ir|  p 


where  c = unit  dyadic  = six’  + yy'  . 


If  a plane  wave  is  normally  incident  on  a circular  aperture  of  radius 


a , then  it  is  evident  that  A = rra.* 


— ^ i and  ^1"  s = Q Th“  —"mmetr*’ 


of  the  circular  aperture  permits  an  arbitrary  polarization  for  the  incident 
wave,  and  it  is  convenient  to  choose  h = y . 

Furthermore,  in  polar  coordinates  the  aperture  will  be  described  by 

Sj  : 0<p<a  0<d<2w 

With  these  specializations,  the  stationary  expression  for  the  transmission 
coefficient  (2)  reduces  to 


t = 


led 


2kira 


a 2s 


a 2m 


A P / f pdpd^zxE  (p)) 

Jo  Jo  Jo  Jo 


a.  2ir 


J'  J p dp  dd  p*  dp  ’ dp' 2xEj^(p).  ,p ').  zxE_j^(p ') 


(4) 


V 
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The  Trial  Function 

Before  choosing  appropriate  trial  functions  for  the  electric  fields  in 
the  aperture,  it  will  be  worth  while  to  examine  some  general  properties 
of  the  aperture  fields.  The  fields  in  question  are  produced  by  a pair  of 
plane  waves  propagating  in  the  directions  2 and  with  magnetic  polar- 
ization vectors  in  the  same  direction  parallel  to  the  y-axis.  The  tangen- 
tial component  of  the  aperture  field  excited  by  the  z-directed  wave  can 
be  assumed  to  be  separable  in  polar  coordinates,  and  to  be  given  by 

£xE^(P,^)  = 2x[pRp(p)ip(^)+aR^(p)|^(^)]  . (5) 

Then,  as  a consequence  of  the  symmetry  of  the  incident  waves  and  of 
the  aperture,  sketched  in  Fig.  2,  it  can  be  seen  that  the  two  tangential 
components  of  the  aperture  field  associated  with  the  z-directed  wave  are 


?.  S_^(p,ir-^)  * E^(p.^)  = Rp(p) 

E_^(p,ir-^)  =4-  E^(p.^)  =-R^(p)  . 


(6) 


The  second  field  can,  therefore,  be  expressed  in  terms  of  the  components 
of  tht  first,  i.e. , 

SxE_^(p,ir-^)  = Sx[^  Rp(p)|p(^)-^R^(p)|^(^)]  . (7) 

If  it  is  assumed  that  the  trial  fields  have  the  same  ^-dependence  as  the  in- 
cident field,  then 

ip(^)  = cos  ^ , f^(^)  = - sin  ^ (8) 

With  (5) , (7)  and  (8)  the  aperture  field  functions  at  the  point  ( p.d)  become 
^ * E^(p,^)  = S X [ p Rp(p)  cos  - % R^(p)  sin  p] 

® ^ * [ ^Rp(p)  cos(w-^)  -!-  Z R^(p)  sin 


(9) 
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so  that 

zxE^(p,^)  = -zxE_^(p,^)  = zxE(p,^).  (10) 

A reasonable  choice  for  the  p -dependence  of  the  trial  field  is  one 

which  satisfies  the  saiixe  boundary  conditions  as  those  satisfied  by  the 

correct  field.  These  conditions  have  been  fxilly  discussed  by  Meixner^^ 
17 

and  Bouwkamp  who  show  that  at  the  rim  of  an  aperture  the  tangential 
component  of  the  electric  field  vanishes  as  and  the  normal  compo- 
nent increases  as  R where  R measures  the  distance  from  the 

field  point  to  the  edge.  Accordingly,  possible  series  representations 
of  the  radial  functions  associated  with  the  aperture  field  components  are 


00 


R (p)  “ ^ 

V a (E.,2n 

/ n'  a ' 

n=0 

oo 

V'  = Vi-(|)^ 

Z -n'l''" 

n=0 

(11) 


where  a_  , b^  are  xmdetermined  coefficients, 
n n 

With  (8)  and  (11),  the  p-  and  ^-components  of  the  trial  functions  for 
the  aperture  field  are  : 


The  corresponding  rectangular  components  are  : 


00 


oo 


^•£(p,^)  = ^ -^1  - ( ^ b^(^^sin^<6 


Vi  (|)^ 
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y-E(p,^)  = [ ^ 


At  the  center  of  the  aperture  these  components  reduce  to 


a^(£;)^+yi  - (|)^  7 a^^l  sin^cos^ 

(13)  ' 


n=0 


$-£(p,l&)  = cos^fS  (a^  + b^)  - 


y-£(p,^)  s (a^ -I- b^ ) sin  ^ cos  ^ 


(14) 


But  since  ^ is  not  defined  at  the  origin,  the  rectangiilar  components 
should  be  independent  of  ^ when  p *=  0 . It  is  necessary,  therefore,  to 
require  that 


a 


o 


(15) 


This  condition  also  insures  tiiat  V-E  = 0 at  ihe  center  of  the  aperture. 

If  the  series  are  truncated  so  that  = 0 for  n^2  , b^  = 0 for  n^l, 

and  the  remaining  coefficients  are  assigned  the  value a^  = -a  , 

a_  = -b^  = 2a  , the  trial  function  reduces  to 
o o 


2 

p E(p,^)  = ■-  P cos  ^ 

(16) 

^•E(p,^)  = - ~ Bin  6 


These  fimctions  constitute  the  low-frequency  exact  solution  obtained 

g 

by  Bouwkamp.  They  are  frequency  independent  and  are  valid  only  for 
ka<^l  . In  reference  (13),  these  functions  are  used  as  a seroth-ordor 
trial  function  by  Levine  and  Schwinger. 

Determination  of  the  Coefficients  in  the  Vector  Trial  Furiction 

The  vector  trial  fimction  z x E(p)  that  was  chosen  in  the  last  section 
contains  tindetermined  frequency-dependent  coefficients.  These  may  be 
chosen  appropriately  by  using  the  stationary  property  of  the  expression 
for  the  tranS^nission  coefficient. 
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If  the  transmission  coefficient  of  (4)  is  written 

t2 

2 


t = — Im  A 
ka^  ^ 


2 im  ^ 


(17) 


ka  4irl . 


then  the  substitution  of  the  explicit  series  for  z x E(p)  shows  at  once  that 


a 2ir  oo 

I,  =r  / pdpd^y  exElp)  (a^B° 

Jq  ^o  n=0 


^ + b B^) 
n n 


r ^ 

= / / P dp  p' dp  ' d<i'  z X E(p).  ; p',^' )•  S X E(p ') 

'^o 


00  00 

= V [a  a + a b + b b ] 

/ /_ n m nm  n m nm  n m nm-* 

n*0  m=0 


where 


a 

>»  . a /■  _ti£ 

" ^oVTW 


2n 


_ h ^ A / n . 2i  *'  ?-i*i 

P v^)  «ip 


n N Py 
'^o 


a'  ' a ‘ 


(£.)2a.£j.j2n 
= / dS  dS'  ^ * 

nm 


■/ 


cos^  cos^'  ® C 


aa 

mn 


ab 

nm 


Cr^  *-/  dS  dS»  {|)^(J^)^8in^cos^T^°5^(p,^;p',^') 

>)  1/T^ 


P^ 


. fis  dS'  (|)^  CO.  t .in*'  P<»> , (p  p '.*')= 

■^1 
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The  are  the  uv  components  of  the  free--space  dyadic  Green's  func- 

tion. With  these  substitutions  the  scattered  amplitude  A^.  takes  the  form: 


A = 


-1 


00 

Z(a  + b B**  Y 

' n n n n ' 


y 4ir  00  00 

Z2 

n=0  m*0 


"V'>'[aa  C^^+ab  C^^+bb 

. ••  n m nm  n m nm  n m nm  ■* 


where  the  B's  and  C's  are,  in  principle,  known  functions  of  the  wave- 
length and  the  size  of  the  aperture,  while  the  a's  and  b's  remain  to  be 
determined. 


dependent  coefficients  and  set 


8A  8A 


Since  the  behaviof  of  the  trial  function  at  the  center  of  the  aperture 

requires  a^  = -b^  extremizing  A with  resnect  to  either  of  these  coef - 
o c-  y 

ficients  gives 

oo 

A >[2a_C*^-(a  -b  )C  -2b 
y X - ni  oir  ' m m'  om  m om-* 

m=0 


oo 


m=0 


(19) 


But  by  definition 
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qj 

= ^I  (a  B“  T b s’*  ) 

2t  1 2iT  / . 'mm  mm' 


m»0 


so  that 


bb  ^ ab 


i i « S’  T»  b» 


[(20““  - C**"  ) a -(2C””-C“")b  ]=r(B“-B“) 
*■ ' om  om ' m ' om  om'  m'  k ' o o' 


m=:0 


A similar  treatment  with  respect  to  a and  b gives 

¥ ^ 


[ 20  ,a  + O^*^  b ] 
*■  vm  m vm  m-* 


= T 

k V 


= 1.2. 


m-u 


^[O***  a + 20**^  b ] B^  H c 1,2, . 

um  m um  m^  K u ^ 


m=0 


Equations  (21),  (22),  and  (23)  constitute  the  linear  set  of  algebraic 
equations  required  to  determine  the  a^  and  b^  . With  these,  the 
transmission  coefficient  can  be  calculated  through  (20). 

In  a practical  csilculation  the  trial  function  must  ha  terminated  at 
some  convenient  value  of  m . The  approadmate  transmission  coeffi- 
cient obtained  in  this  way  will  be 


tW)  . .. 


»m»ml  - 


msO 


m«0 


where  N = j + k is  the  order  of  the  approximation.  For  a first-order 
solution  with  j * 1 , k ■ 0 , the  set  of  equations  to  be  solved  consists 
only  of  (21)  and  (22)  with  vs  1 : 


S + c‘j>)  aj  4(B»  - B^)  (25) 


5 
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I 

i 


These  equations  are  readily  solved  for  the  coefficients  a^  and  , 
which  can  be  substituted^ into  (24).  The  result  is 


t ^ iKi  - ^Qta)  + R(v^) 

^ Q^(a}  - 4P(o)R(a) 


(27) 


Here,  a = ka  and 


P(a)  = C“ 


Q(a)  = 2C5“  C„^J> 


R(a)  . C“  - C*o‘>  + Co‘-*>  . 

These  integrals  are  evaluated  in  Appendix  A.  In  additicn, 

a 


a 

'o 


B-  = waT  -£-S£_ 


= ir  a 


a 

'1  a 


J \ri — r 

o y A - p 


2 ,2 
IT  a 


- r — T 

a p 


/ p/ 


ip  - - 


1 _ 2 

« ci 


The  numerical  result  of  t^^^  is  compared  with  previous  results  in 
Fig.  3.  It  can  be  seen  that  the  variational  approximations  are  quite  close 
to  the  exact  solution  and  that  t^^^  provides  a much  higher  degree  of 
accuracy  than  in  the  range  1.  5^ka^4.  5 . 

A further  comparison  of  interest  can  be  made  for  small  values  of 
ka.  If  the  transmission  coefficients  obtained  by  the  various  theories  are 
expanded  into  power  series  in  ka  , the  re  stilts  are  : 

g 

(i)  £xact  Solution  (due  to  Bouwkamp  ) 

4 

tg  = [ 1 +22  ^ Q 39793197  (ka)**  + . . . ] (28) 

27ir^ 


1.8 
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(ii)  Zerotlx-order  approximation  by  the  variational  method  (due  to 

13 

Levine  and  Schwinger) 


(iii) 


27ir 

First-order  s^proximation  by  the  variational  method  obtained 
in  this  report 


[ 1 + II  (ka)^  + 0. 40790023  (ka)'^  + . . . ] (29) 


t(i)  = 1 1 + II  (ka)^  + 0.  39678912  (ka)^  + . . ] (30) 


27x 


tt0>. 


The  third  terms  in  the  bracket  show  that  t'  ' is  more  accurate  than 


U) 


III. 


The  Elliptical  Aperture 


If  the  aperture  in  the  plane  screen  is  elliptical  in  shape,  the  general 
procedure  of  section  11  is  much  more  difficult  to  apply.  Either  elliptic  or 
oblate  spheroidal  coordinates  must  be  used,  and  the  computational  labor 
would  be  ^eatly  increased.  However,  considerable  improvement  in  the 
Kirchhoff  result  can  be  effected  by  using  a single-component  trial  function 
for  the  electric  field  in  the  aperture  corresponding  to  the  original  trial 


For  a normally  incident  plane  wave  polarized  so  that  E^°^^(r)  =$e 


ikz 


<uia  ^ - 


jjinC(^j  _ ^ ^ _£  fvinction  is 

1 

E(p)  = ^(l-p^)"^ 


(31) 


Substitution  of  this  function  into  the  stationary  expression  for  the 
transmission  coefficient  (4)  gives 


-2ir  ab 

■“lir- 


Im 


(32) 


1 2x 


= r r p do  p'  dp'  d^'Yl  - Vl  - p'^ 
o o 


where 
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cx> 


( p P 


dk  dk 

X 


^ (1  - erp  [ikx(x-x')  + ikyijr-y')], 

^ 1/1_2  i_  i 2 ^ 

-CO  1/k  -k  -k 


lA 


X V 


the  integral  can  be  simplified  to 


I ■ 


4k‘ 


2f  oo 

II 


dv  doV^^^l  j|/,  (kLv) 


O ''O 


T 3 2 

Li  V 


2f  CO 

// 


r2 


dv  dO 


cos  0 


o o 


3yCI 


L^Vv“-  1 


(34) 


where 


L = V (a  coo  0 )^  + ( b sin  0 


The  procedure  is  similar  to  that  developed  in  reference  (14).  Finally,  by 
carrying  out  the  integration  with  respect  to  v , we  find  that  e 


I = 


7 

a o H / 
/ 


4k‘ 


[F(kL)  - coo^O  G(kL)]  ^ 


(35) 


where 


F(x) 


1 

-G 


2x 


2x 


1 + 
w 


S^(2x) + Sj(2x)  - (1  + ^) 


/ 


S„(t)( 


*The  first  integral  is  related  to  F.  ,(kL)  in  reference  (19);  the  second  is 
related  to  l,(kL)  in  appendix  A. 
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G(x)  = 


J 

x|^  2x  ' 


(2x)  - Jj(2x)  + (-2  + -^ 

4x 


ex 

r 


J 


J„(t)  dt] 


2x 

+ S„(2*)-S,(2x)  + (4*^>  J S„(t)dt]| 


Numerical  values  of  the  transmission  coefficient  are  given  in  Fig.  4 
for  0<ka<10,  and  b/a  = 1/2,  1/3.  These  are  plotted  together  with  the 
corresponding  result  for  the  circular  aperture  which  has  been  given  in 
reference  (13),  and  which  can  be  obtained  from  equation  (32)  by  pvitting 
b/a  = 1 . Using  the  e3q>ansions 


1 2 ^ 
1 X . X 


F(x)  s -^  + ^+  . . .]  + 


1 2,4 

1.x  3x 


4x^  16x^  . 32x^ 


.] 

(36) 

. .] 


in  (35),  a form  of  the  transmission  coefficient  appropriate  to  small  values 
of  ka  can  be  obtained  : 

t = (J)  (ka)^  [A.(^)  + A,(^)  (ka>^  + A.(-)  (ka)^f  . . .] 


^Ytr 

(37) 

b b D 

where  A (— ) , A, ( — ) , A^(° ) ....  are  functions  of  the  eccentricity. 

Q sL  ^ A A 

They  are  given  in  appendix  B.  With  small  eccentricity  (37)  becomes 

♦ , 64  , a . ^.4  r , , 9 Z v . / 27  349  2 . i yi,  v- 

t = 2(^)(ka)  l(l-^e  . . *)'*'( ~ 2TSI!F  c +•  • • ) (^) 

27it 


. y 8937  155223  2 . , „ v4  . , 

'TIZW  " 4^00  ® +•  • •)  (Ita)  +.  . . J 


(38) 
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Appendix  A 

In  this  appendix  we  will  sketch  the  method  used  to  evaluate  the  ftinc  - 
tions  required  for  the  transmission  coefficient  of  the  circular  aperture. 
These  functions  are  given  hy  (27): 

P(tt)  = 

Q(q)  = 2CQ^^(a)  - 

R (a)  = C “(a)  - 


The  integrals  czui  be  simplified  by  making  use  of  an  integral 

representation  and  an  addition  theorem;  the  integral  representation  is 


20 


00 


^ik  j / u)dfa) 


where  arg  (to^  - k^ ) ^ 


V 

"T 


ta)>k 

(i)<;k 


and  R = (p^  + p’^  - 2 pp'cos(^  - ; 


iumziOu  tneOfein  is 


(2  - 6^)  Jj^(cop)  J^(wp')  cos[n(^  - , 

n=0 


«on=  0 


for  11  = 0 
for  n ^ 0 


With  these  substitutions  (18)  gives 


TR164 


-17- 


1 t (I)"” 

^rnn  ' 3?  I 


aa 


4>  y 1 -{^)‘ 


a ^ p * ^2n  2ir  2tr 

r ■ ^ p'  dp'  C 008  6 f cos  ^ 

o y 1 -{^)^  -^o  4) 


oo 


[C03{^-^)--^  7 I?  Spn)— 


n=0 


oo 

f 


'o  ,VV-k 

The  orthogonal  properties  of  the  trigonometric  functions  permit  the 
integration  with  respect  to  ^ and  to  be  carried  out.  The  result  is 


. aa 
'nm 


X Vu^-k^  X^i  ■Joy 


» ,41,2m 


^ 


p'  dp’  [ J (wp)  J («p') 


o^oyi-(i^) 


+ J2(up)  Jj,  - “2 


2 1 


PP 


(up)  Jj  (up')] 


By  introducing  the  notation 


= I 

on 


r (|;)^  P dp 


Jo  yi-(|)' 


Jq  (up)  . 


a r ^a^^  P ^P 

l2n  = / 

•4  l/i-(i 


yi-(|)' 


J2  («p) 


the  integral  can  be  written  as 


. aa 
'nm 


00 

ir  r u 


du 


k‘ 


a .a 

om 


+ 1^1^  ] +i4[I^  -if  ][I  * -I  * 
2n  2m-*  j^2  . on  2n^  '■  om  2irj 


•') 
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01 


a 3t  / * j ir/2 

p JJwa)  dp  r 2 3 

a sin  0 J^(uasinO)  dO 


/ kfl T 

>o  a ya  - p 


■/ 


-a  T ^ ^ 

■ ^00  ■ ^00  ■ ^ 


Vl  ^T~yrJ  - T~  W^ 


(«a) 


(«a) 


'20 


a 


a p J dp 

'-'o  ]fL~^  - p 


7 T 


viz 

= j a^  sin  0 (wa  sin  0) 


dO 


viz 

2 ^ 2 

= a / f»inO  [ — ^ J,  ((uasinO)  “ J (ua  sin  0)1  dO 

/ {oa  sin^  1 ' 'o' 


oja  sin  I 


2^/ir  r ^^2^“^^  \lZ 

ITT 


’ “T?  ‘ ~.’5/J 

V“*“i 


(wa) 


'20 


a l/72  Z”  - I » . ir/2 

r P ya  -p  Jylup)  ap  r 2 2 

= I / a cos  Osin 0 J2(uasinO)dO 


= I?.  - :f.  = a^1/J  r 2 - 


V,,  ■W'"’ 


£*1 


r 6 


(wa) 


- 3 


■^V2<”'  , 
(«a)  ' 


These  forms  for  the  I*e  are  obtained  by  maMng  use  of  Sonine's  first 
finite  integral 


T ^ ^ ST 

“n+v+1 ' ^ 


2"r(v+i) 


Re  |i  ,v  > -1 


and  its  alternate  form 


J (z  sinO)  sin 
n'  ' 


1- 


‘‘O 


cos 


2V+1 


0 dO  = 


° ti-f  V,  V 


. Rev  >-l 
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where  s (z)  is  a Lonunel  function, 
these  functions  are  defined  by 


When  the  indices  are  half  integers 


Sji  „ (z)  = ^ sin|iii[  J^(z)  / z**  J_^(z)  dz  - J_^(z)  J^(z)  dz] 


and  are  related  to  the  Struve  function  S^(z)  in  accordance  with  the  equation 

V, Hi.  s^(«) 


22 


When  the  integrated  functions , 1 , are  substituted  into  the  expressions  for  the 
C's  and  these  are  combined  according  to  equation  (27),  a certain  amount  of 
algebraic  manipulation  yields 


P(o) 


, . 00 
v a / {•> 

IF. 


db) 


(b>a)  (coa)"‘  (wa)"^  (bja)"* 


Qio)  = 


P . oo 

= ^ / r 

-L  V u^-k‘ 


, Jl/z(“»>  , V--®’  ■'V2‘“>  . „ . J 

2 — -b  “ — > +H » 4—2  5 — 

(«a)  («a)^  (oa)'*  k""  (wa)^ 


o - °P 


R(a)  = 


tr''a'  I fa)  dfa> 


'■'o  Vfa)  -k 


T 2 / 

~l/Z'~ — ' 
(«a) 


T I . V T / . ...  \ 

-2-=i-= + 5 


(«a) 


t2 

(«a)^ 


> T 2 

"3/2^  • 

k \wck) 


It  can  be  seen  that  each  of  these  functions  is  a linear  combination  of  the 
following  integrals 

, _ T .i.  4<“»> 

‘ -i 

I f ^3/2<“> 

' Jo  k"i/7:j2 
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oo 


I.  = 


Jo 


ci>  d(t)  . , <•>  \ 

11 T I 


(•a)^ 


oo 


iii  dci) 


/7TT2- 


Va)‘'-k 


J^/^(u>a)  J3y^(u)a) 

(«a)^ 


The  remainder  of  this  appendix  will  be  concerned  with  evaluation  cl 
these  integrals. 

The  integral  1^  can  be  rewritten  as 


I = i f ^ 1 

1 a 1 . /— 1« — ^ a 


a“) 


/ V 2 i 
'a  r s*>  -a 


By  making  use  of  the  integral  representation 


Jo  y 

23 


2 Z" 

a -w 


J^(z)J^(z)  = 


ir/2 

'/O 


J . (2z  cos  0)  cos  (u  "V  ) 0 dO 

Jl  + V 


Re  (^.+v  )>  -1 


The  real  part  of  becomes 

it/2  00 


j r f J,(2cjcos  0)  doi 
Re  Ij  = ^ / d«  J 


But  with  the  infinite  integral  formula 


J \rr-T 

/ o>  -a 
24 


J («a)  d'u> 


v-1 


a>0 


Re(^-^)>Re  n>- 


the  integration  with  respect  to  u can  be  carried  out  and  the  result  is 
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T 12 

1 ^ 


J,,_  (2a  cosO) 


y ia  cos  d 


We  now  make  use  of  the  fact  that  when  v is  half  of  an  odd  integer,  the 
function  J (s)  has  a finite  representation  in  terms  of  algebraic  and  trigo- 


25 

metric  fxinctions  of  z , i.e. 


&f  V 


-V**  “ * 


ir/2 


Then 


J. 


^ i F,(a)  . 


tra  a 


where 


v/Z 


Fj(o)  = 


3m  cos  9)  do 
Z COB  6 


It  is  evident  that  Fj^(O)  = 0 and  that 

ir/2 


*^(o)  = j cos  (2a  cos  0)  dO  = ^ J^(2a) 


From  these  it  follows  that 


Fi(a) 


ca 

F'(a)  da  = 5 / Jo(2«)  = J / 3 J,t)  dt  . 

Jo 


'J 


Accordingly, 


Re  I,  = ^ / J (t)  dt  . 
1 2aa  / o'  ' 


The  imaginary  part  of  can  be  expressed  by 
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Im  I, 


f ± 7 a« 


(Zojcos  0)  dw 


o K a -Cl) 


By  making  use  of  the  second  form  of  the  Sonine's  first  finite  integral^  , 
the  integration  with  respect  to  co  can  be  carried  out: 


2 r ®i/2  1/7(20  cos  0) 

Iml,  = -=-  / dO  

y 2o  cos  0 


Therefore  the  single  integral  can  be  rewritten  as 


ir/Z 

■‘■=2 


[1  •cos(2acos0)]  dO  2 \ 

^ *^2'*^' 


where 


= / 


1 - cos  (2a  cos  0)1  d9 
2 cos  0 


Again,  with  F,(0)  = 0 and  with  the  integra.1  representation  for  the 

26  ** 

Struve  function 


S (z) 

V 


2(7)'' 


r(v  + 7)p(4) 


sin  (z  cos  0)  sin^^O  dO 


we  have 


Re  V >- 


a 

r 

F'^Co)  = I sin(ikico8  0)  dO 


= 7 S^(2o) 
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Cons  equently, 


ra  2a 

(a)  dn  = ^ I'  (2a)  da  = J T S^(t)  dt. 

ft  '/o 


and 


2a 

^ ^ laS  I • 


Therefore 


~ 2aa 


2o 

1 r r 


2a 

r 


[ I j^(t)  dt  + i / s^(t)  dt]  . 

Uq  Jo 


The  integral  1^  can  be  evaluated  by  decomposing  it  into  two  parts: 


00  ^ 

Hi 


dw  •'3/2 


(«a) 


dwJ3^2^“> 


1 r ^“•'V; 

J 

'Jo  r w -< 


If  we  let 


oo 

r 


do)J3^2t«) 


aa  J V ”Z  Z" 
‘'o  r « -a 


- ^-1/2<“>1 


J _ r dh.j^^(u>) 

® ~ / if  n 

'^o  (i)f  cj  -a 


'6' 


00 

•i 


d«  Jy2^w) 

V 2“Z’ 

y w -a 


then 


^2  = 

aa 
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The  evaluation  of  Ig  proceeds  in  a rnannitr  sinciilar  to  the  one  used  for 


Ij.  Thus 


^ J-y,((u)  J,y,(w)  d«  , ^ J,(2uco8  0)  dw 

Re  = 4 / cos  0 do  / -= 

' \rr~T  ^ I J \rr-T 

f Ir  « -a  u)  Kw  -a 


a u 


ir/2 

= 1/  co»0d0j/| 


^cos  6 d 


T ^3/2(20  008  6) 


T(lZ 

" f ' 

-A 


- CO.  (2.  «.<.)) 


and , 


2a 

= [ i Fj(o)  - 2 Jo(2^*>]  " f " 2J^(2a)] 

ira  4a 


Im  I, 


^ ^ JyzCo)  Ji/2(.o)do.  2 ^Jz<^co.g) 


dci) 


•nIZ 

c.  i cos  0 dw  , ^ 


But 


1/2  2 

*3/2  -3/2^*^  = z [ z - 2 sin  z + — (1  - cos  z)] 


Therefore, 


Im  Ig  = 


w/2 

o/„ 


dS  faco.«  - .in(2aco.g)  + ii--.  co.  <^co.  gj],  , 

' ' 2acos  6 ^ 
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Similarly 

Re 

where 

F^(a) 
Since  F»(0) 

F^(a) 

therefore 

F3(a) 

Re 


tra 


— S'  [— 0-2S  (2a)  + 
. 2 '■  IT  o' 

4a 


2a 

- W 


S^(t)dt] 


tt> 


■1 


i/  t i 

y oj  - a 


ir/2 


j jr  ^ J.(2wcos  0) dcj 

= - / cos  20  dO  / — 

Vo  j.  fJTT’ 


viz 


-§f  coBZOdSyi 


y 2a  cos  8 


"1/z  (2acos  Q) 


ir/2 


- 2 f 

’ / 

'''O 


20  dO 

2a  cos  0 


= ^ *-3 ‘“I 


ir/2 

r 


= 0 and 


■tr/2 


cos  20  cos  (2a  cos  0)  dO  = ~ ^ 


a 2a  2o 

(a)  do  = - 5 r ijd)  dt  = - 5 r j^(t)  dt 

'Jn 


2a 

'-'n 


= it  + I <2a) 


and 


M=* 


s 
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On  the  other  hand 


^ V ~7f=i  ^ ” 

/ a -w 


cos  20  dO 


Jj^  (Zujcos  0)  dO 

irz~T 

r a -0) 


= - / cos  20  [ ^ - C08(2acos0)]  dO 

^ 2a  cos  0 


2 a COS  0 


where  F^(a) 


= sin^e  I ^ - cos  (2a  cos  0)]  dO 
J 2 cos  0 


Again  (0)  = 0 


IT/ 6 

F^  (»)  = sin^O  sin  (2o  cos  0)  dO  = S^(2a)  ; 


therefore 


F4(a)  = ^ / 


_ ^ S,  (t)  dt 


If  the  recurrence  formula  for  the  Struve  function  is  used 


® S'y(z)  = - vSy{z)  + z S^_j(z) 


we  have 


2o 

^4(0)  = \{  f S^(t)dt  - Sj(2a)] 

'•/ft 


Consequently  Im  1, 


<ca 

= 1 2Sl(2“)  - f S„(‘)  dt  ] 

>Jes 


-Ml#’ 


'll 
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Combining.  1^  and  1^  , we  obtain 

^2  = :b-  '*5  - U 


ab 


J. 

aa 


2a 

1 I * f y*-*'  - I ^ V 


2a)] 


The  next  integral  to  be  considered  is 


‘■■7 


u>  db)  I , . *^3/2 

— (i“7z)  


^ V 2 , i k“  , ,3 

o It  fa)  -k  (b)a) 


2 — F 2 

“ J^2  (<-a) 


With  an  integration  by  parts 


r iTTTr 


But 


therefore 


“ J^/,(u>a) 


aa 


i2a/  Jy  (coa)  (oa)  da.  + / -- 

-'o  " yiTik^ 


3 *^3/2^ 

^/2<»>  = - I " Jl/2<">  • 


oc 

[3  f 

Jo 


00 


dfa)  ^ Vtd^-k*'  Jj^2^wa)dfa> 


- 2a 


fa) 


J. 


(t) 
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and 


■7 


3^^(ua)  d<i) 


] = 


o Yw^-k 


I-  = 


- -I  ^7  - ^2 
aa 


2aa 


. T . 1 T 

2 "7  ? ^2 


where 


^7 


a ^ y w^-a^  -^3/2^“^ 

7 

Jo 


But  it  is  easy  to  show  that 


2a 

= ~ f 

4a  / 


yt)dt+i[-^  + ;^  i s^(t)dt! 


Therefore  when  and  ly  are  combined,  we  find  that 


-1 


aa 


2a 

r 


4o  / 8a 


2a 

4a  / 8a 


+ * 


The  last  integral,  , is  directly  related  to  1^  as  follows 


00 

■I 


udbi  = - r ^ j 

'/o 


(<«) 


a o 


Therefore 


“4 


4aa' 


[ 


ij 


J^(t)dt-  2y2a)]+i 


2a 


S„(t)dt  - 2S^(2a)] 


\ 
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Finally,  by  subatituiing  the  integrals  Ij  , into  the  fvinction 

P (a)  Q (a)  and  R (a)  , we  obtain 


2 3 
ir  a 


P(a)  = ^ 
4a 

2a 

I 


r f •%  <XX  X X X.  XX  X k 

^(2^ -g- _j)  Jo(2q)  ~ - Jj{2a)  +^--^_  + 

a a 


J„(t)dtl  “ i + ,2-  “ J^)  S„<2a)-^iSj(2a> 

a 


+ ( 


2a 

a 1 1 , 11  1 , r (t)  dt] 

1-  f -i  * TS  -%>  J ° 

A 


Q(«)  = 


[(3--ij)  J^(2o)  - I Jj(2a)  + (o 


2a 


j^(t)  dt] 


2a 


- I J * / S„  (t)  dt] 


R(a)  = 


ir^a’ 

XT' 


2a 


® “ o/„ 


a 


_ i 


^<?-u 


2a 


Sj,(t)  dt] 
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Appendix  B 

In  this  appendix,  a form  of  the  transmission  coefficient  for  the  ellipti- 
cal apertures  which  is  appropriate  to  small  values  of  ka  will  be  derived. 

When  the  expansions  (36)  of  the  functions  F(x)  and  G(x)  are  substi- 
tuted into  (35),  and  the  i!itegr:»iions  with  respect  to  0 are  carried  out  term 
by  term,  we  find  that 


2 

I = [R  (e)  + R (e)  (ka)^  + R,(e)  (ka)^  + • • • 

4k^a  « ^ 4 


+ i [ye)  + 12(e)  (ka)“  + ye)  (ka)'^  + • • • ] 


where  B.g(e)  = [ Sl-e)  ••  K {e)j 


= ire  K(e)  -i-7  [K(e)-E(e)] 
15e 


^4^*^  ~ “IT  '*■  ^ > [E(e)  -?(e)] 

35e 


I,(e)  = 


- d r o 2-t 

[3  - e 3 


^4^®^  ~ ^5125  [ 16  - 10  e^  + 3 e^] 


e = eccentricity  = y 1 - (-j)* 

K(e)  = complete  elliptic  integral  of  the  first  kind^® 


ir/2 

I ^ 

l/l  - e^ 


sin^O 
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_ wr.  .^Iv2.2  . f l 3v2  4 . .1-3-5.2  6 , 

e 4(^)  e ^ijrris)  c + 'J 


E(e)  = complete  elliptic  integral  of  the  second  kind 


ir/Z 


1 - e sin  0 d9 
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Substitution  of  I into  (32)  gives 
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where  a = ka 
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where  A*(e)  = 
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